Geometric phase with nonunitary evolution in presence of a quantum critical bath 
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Geometric phases, arising from cyclic evolutions in a curved parameter space, appear in a wealth 
of physical settings. Recently, and largely motivated by the need of an experimentally realistic 
definition for quantum computing applications, the quantum geometric phase was generalized to 
open systems. The definition takes a kinematical approach, with an initial state that is evolved 
cyclically but coupled to an environment — leading to a correction of the geometric phase with 
respect to the uncoupled case. We obtain this correction by measuring the nonunitary evolution 
of the reduced density matrix of a spin one-half coupled to an environment. In particular, we 
consider a bath that can be tuned near a quantum phase transition, and demonstrate how the 
criticality information imprinted in the decoherence factor translates into the geometric phase. The 
experiments are done with a NMR quantum simulator, in which the critical environment is modeled 
using a one-qubit system. 



PACS numbers: 

For decades, the geometric phase [I] (GP) has fasci- 
nated physicists for its elegant theoretical grounds and 
its practical applications . The GP is resilient to dy- 
namical perturbations, thus, it might serve as a natu- 
rally fault-tolerant quantum information processing de- 
vice 3J. In order to explore such applications, and un- 
like traditional studies of the GP in closed systems with 
pure states, one must take into account realistic exper- 
imental conditions — i.e. the explicit presence of noise 
and environments. Uhlmann was the first in consider- 
ing a system in a mixed state, embeded, as a subsystem, 
in a larger system that is in a pure state [Ij. Later, 
Sjoqvist et al [5] put forward a definition of the GP for 
a general mixed state undergoing a cyclic unitary evo- 
lution — subsequently measured using NMR interferom- 
etry in Ref. [6]. Different approaches to the problem 
were proposed [7j. In the present Letter, we will follow 
the line of Tong et al. [Q, who developed a kinematic 
generalization of the GP to open systems that takes into 
account the coupling to an environment (leading to a 
nonunitary evolution of the reduced density matrix of 
the system [5].) Arguably, this approach is better suited 
to explore the usefulness of the GP in a real quantum 
computer undergoing decoherence processes ^0\. Here 
we report a measurement of the GP for a spin 1/2 un- 
dergoing nonunitary evolution induced by the coupling to 
an environment, using the decoherence factor or fidelity 
decay [llj . In particular — motivated by the recent obser- 
vation that baths near quantum criticality induce strong 
decoherence [llj — we choose an environment that can 
be tuned near a quantum phase transition (QPT). This 
choice not only adds richness to the behavior of the GP, 
but also advances the program of understanding it in gen- 
eral open systems jlOj . In our experiments, performed in 



a NMR quantum simulator, we measure the full time de- 
pendence of the decoherence factor of the system-spin 
— from which we can determine the GP using the re- 
sults of Ref. For the environment, we introduce a 
simple qualitative model of the ground state degeneracy 
that occurs at QPTs. Apart from demonstrating an al- 
ternative to traditional interferometry-based approaches 
for measuring the GP in open systems, our results fur- 
ther establish the strong connections between quantum 
information, quantum criticality, decoherence, and the 
quantum geometric phase [TTI - fT^ that have been the fo- 
cus of much recent research (especially in the context of 
quantum simulations [El [16]). 

The correction to the GP by a critical environment was 
first studied by Yi and Wang [T7], who gave some gen- 
eral analytical results and found numerical instabilities 
in the GP of a qubit near criticallity of the bath (an XY 
spin chain). More recently, it was shown that the GP 
of a spin coupled to an antiferromagnetic environment 
changes suddenly when the bath undergoes a first order 
QPT [18j. Notice that our problem is seemingly related 
to, but different than, the use of the GP as an order pa- 
rameter in a QPT of a closed system, as studied first by 
CaroUo and Pachos and others [T2l [T4] . 

We consider a spin 1 /2 coupled to an environment with 
a total Hamiltonian H = flZg ®l£ + Zs® Hss +Is®H£, 
where Hg is the Hamiltonian of the bath, Zg is the 
z Pauli matrix of the system. Is is the identity oper- 
ator of the system and Ig the one of the bath. For 
simplicity, we only consider a dephasing spin-bath in- 
teraction, Zs ® Hs£, neglecting relaxation effects and 
limiting the relevance of the initial state (see discus- 
sion below). We take a product initial state for the 
spin-bath system, p(0) = It/iq) (V-'oI ® |£(0)) (£(0)|, where 
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IV'o) = sm{0/2) |0) + cos(6'/2) |1) and |e(0)) is a general 
initial state of the bath. In absence of the bath, the spin 
follows an evolution around the Bloch sphere, reaching 
again the initial state for t = 2tt /VI. To compute the 
global phase gain during the evolution, one can use the 
Pancharatnam's definition [19 , which contains a gauge 
dependent part (i.e a dynamical phase $^ = — 7rcos(0)) 
and a gauge independent part, commonly known as geo- 
metric phase $g = 7r(l — cos{6)). When coupled to the 
bath, the reduced density matrix of the system at a time 
t is 

p^{t) = sin2(0/2) |0) (0| + cos2(0/2) |1) (1| 

+ '^e-^'^'r(t) |0> (1| + ^e^^^V* W |1) (0| . (1) 



Here, r{t) — (eo(i)|ei(0) the decoherence fac- 
tor induced by the environment, with \ek{t)) = 
^-it[HB+{-i)''HsE] \s[Q)). The phase $ acquired by the 
open system after a period r is defined as |5], 



(a) 



arg 



(2) 

where \k{t)) and efe(T) are respectively the instantaneous 
eigenvectors and eigenvalues of Pi-(t). Of the two k modes 
(+ and — ) of the one qubit model we are treating, only 
the + mode contributes to the GP. Because our envi- 
ronments can induce a complex decoherence factor, i.e. 
r(t) = |r(t)|e~*'^, we obtain a slightly different expression 
than that shown in Ref. 9], namely 



tan 



cos(0+/2) 
sm{0+/2) 



(3) 

(4) 
(5) 



sin(^(r) sin(^) sin(|) 
cosi^(t) sin(^) sin(|) -I- cos(^) cos(|) 
where we have defined 

2(e+ -sin2(6'/2)) 
y'|r(t)|2 sin^e) + 4(e+ - sm^{9/2))^ 

\r{t)\sm{9) 
^|r(t)|2sin2(0)+4(e+-sin2(0/2))2 

During normal quantum evolution, the system gains a 
global phase. The central result of Eq. ([s]) is to extract 
(by a proper choice of the "parallel transport condition" ) 
the purification independent part of the phase — which 
can be termed a geometric phase because it is gauge in- 
variant and reduces to the known results in the limit of 
a unitary evolution. 

We have studied Eq. ^ both numerically and ana- 
lytically with an Ising spin chain environment (see sup- 
plementary material [21]), a paradigmatic example of 
a quantum phase transition. In particular, the non- 
analiticity of the GP at the critical point in the ther- 
modynamical limit becomes evident in the limit of weak 
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FIG. 1: (a) Observed absolute value squared of the deco- 
herence factor and (b) its argument, both as a function of 
time and external magnetic field strength, (c) Computed 
correction to the geometric phase for a choice of 9 — n /4. In 
large filled circles is the experimental data, and in solid line 
is the theoretically expected value (without free parameters). 
The corrected GP is the difference between the GP measured 
in presence of the environment (small filled circles), and the 
GP measured when the system and environment are decoupled 
(small empty circles). In this setup, = IOOtt Hz, A — 0.0251, 
and S = O.lfi. 



system-bath coupling |2T]. Nevertheless, a full quantum 
simulation of a large enough critical system is on the 
edge of current technology, and beyond the scope of this 
Letter. Therefore, before turning to the experimental 
results, let us discuss briefly our choice for modeling a 
critical environment. 

Near its critical point, the spectrum of a quantum crit- 
ical system is characterized by the closing of the energy 
gap between the ground and the first excited state. In the 
thermodynamical limit, the gap closes with a power law 
^ |A — Acl^'' (where z and are critical exponents), but 
for all finite size systems there is a minimum gap A near 
Ac- It is remarkable that, for many purposes, this feature 
of the spectrum is enough to describe qualitatively the ef- 
fects of a critical environment: as long as the excitations 
involved are small, and one is only interested in qualita- 
tive behavior, a small energy expansion of the decoher- 
ence factor can justify considering just two levels with 
appropriate dynamics |16j . Thus, we propose to mimic 
a complex critical bath using a simple two-level system 
model with Hamiltonian Hg = A|A|^'^ ^ A + A Xg, 
where Ac = represents the "critical point" or QPT. 
The simplification might seem excessive, but it has been 
used successfully before: For ziy = 1 (the mean field ex- 
ponents), it gives a correct qualitative description of the 
creation of topological excitations during a finite speed 
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quench [52]. In essence, the model is quantitatively not 
far away from the small systems used in demonstrations 
of quantum phase transitions [151 [53] . A complete char- 
acterization of when this model does not describe the 
correct physics of a critical model is missing (one such 
example is the calculation of the path length of an adia- 
batic evolution [H]). Nevertheless, our results show that 
for the GP problem the model gives a fair description 
when the gap A is much smaller than the natural fre- 
quency fl of the system spin. 

Using a tomographic approach, we measure the GP 
of a qubit coupled to a critical environment using a 
nuclear magnetic resonance (NMR) quantum simula- 
tor, with the environment represented by the two level 
model described above (with critical exponents zv = 1 
and a dimensionless transverse field strength B = AA). 
The target Hamiltonian to simulate experimentally is: 
H = nZs + SZsZe + BZs + AXg, where Zs and Zg 
are the z Pauli matrices of the system and environment 
respectively. We obtain the GP by measuring the mag- 
netization of the system spin in the X — Y plane, which 
gives us the decoherence factor r{t). 

Denoting by e± = iAVlTA^^ the eigenenergies of 
Hg, the decoherence factor of this model is 

r{t) = e"-(^)* [cose_(A + (5)<- 
' ^r(A)MAW""-^^n- 

where, to simplify notation, we have chosen the system- 
environment interaction to be Hs£ = S{Is — Zs)Z£. 
The correction to the GP due to this decoherence fac- 
tor [shown in Fig. ([I]) with the experimental results to be 
discussed below] contains the main elements observed in 
more complex models, as Ising spin chains [21] : a max- 
imum correction of the GP at criticality, and a small 
asymmetric correction far away from the critical point. 
From this simplified model we can get insight into the 
physics of true quantum critical baths. 

Overall, the experimental sequence is as follows: We 
first fix ri, 6, and A. Then, for each value of B, we 
initialize the system, and measure the decoherence fac- 
tor r{t) of the system after evolution with an operator 
U = e~*'^* for various times t £ [0, 27r/17]. The measured 
decoherence factor is shown in Fig. and ([!]}). The 
GP is calculated using a numerical interpolation of r{t) 
in Eq. ^. 

We choose the C^"^ and spins in the molecule of 
carbon-13 labelled chloroform (CHCI3) dissolved in d6- 
acetone as the quantum registers (qubits) for the demon- 
stration. The C^'^ atom simulates the system, and the 
the environment, where the scalar coupling between 
them is measured to be J = 215 Hz. Data were taken 
with a Bruker DRX 700 MHz spectrometer. 

Our choice of system-environment coupling makes the 
decoherence factor r{t) independent of the initial state 
of the system (given by the angle 9) [see Eq. (1)]. This, 
in turn, makes the GP depend trivially on 6, which can 



be fully appreciated when approximating Eq. (|3| in the 
weak coupling regime (see supplementary material [2Tj). 
Because we concentrate on how the criticality properties 
of the bath affect the GP, it is experimentally convenient 
to fix an initial state of the system that maximizes the 
signal to noise ratio, and change only the parameters of 
the environment spin. In particular, we choose the input 
state of the system to be ([0)5 -I- \l)s)/V2- The corre- 
sponding decoherence factor can then be used to compute 
Eq. ((3| for any other initial state of the system. From 
Eq. (Tj we can see that r{t) is encoded in the coherent 
terms proportional to (2(t+) [see Fig. (26)], which can 
be observed directly in NMR by adding the two complex 
amplitudes of the peaks in the C^^ spectra. 

We use the gate sequence of Fig. ([2]-a) [55H27] to 
prepare the pseudo-pure state |00)5£, to which we ap- 
ply the unitary g-^'^^s/4:giaY£/2 reach the input state 
\An) = (|0) + \l))s\g)£/V2. Here \g)£ is the ground 
state of the environment for a given S-value, \g)£ — 
|0)cos(q!/2) - |l)sin(a/2), where tana = -A/B. Be- 
cause the decoherence factor is independent of the initial 
state of the system, we chose it such that it maximizes 
the signal-to-noise ratio of the experiment. 

The quantum simulated evolution U for a time t can 
be implemented with a Trotter approximation [28l [29] , 

jj ^ ^-iAtXe /2^~iStZsZe g-intZs ^-iBtZe ^-lAtXe /2 ^j-^ 

where we choose il = IOOtt Hz, A — 0.02il, 5 — 
0.117, and we apply the evolution operator r/t times. 
We checked numerically that the Trotter approxima- 
tion reduces the fidelity less than 0.3% for B e 
[— 0.217, 0.2ri]. Furthermore, we decompose the uni- 
tary operations e"'^*^^ as e-'^^^^'^e-'^^^^e''^^^/^, and 
^-zQtZs as g-^-^Xs/4^-tntYs^i7rXs/4: ^^^^^ imple- 
ment them with standard rf pulses. The coupling opera- 
tion e-"5*-Zs2£ jg realized using the natural spin coupling 
with an evolution time 25t/{'KJ). After the evolution J7, 
we measure the magnetization in the XY plane, which is 
proportional to the decoherence factor r{t). The whole 
gate sequence for each measurement is shown in Fig. ^ 
b). Notice that we measure the absolute value as well 
as the complex phase of r(t), necessary for the GP. The 
total evolution time was always well below the natural 
decoherence time of the quantum simulator. 

To eliminate systematic errors, we repeat the exper- 
iment but uncoupling the system and the environment 
((5 = 0). From this we compute a baseline GP, which we 
subtract from the full (coupled) experiment. Thus, we 
obtain the correction to the GP due to the presence of 
the critical environment, which agrees well with theoret- 
ical expectations [see Fig. Qc]. 

Conclusions. Using a NMR quantum simulator, we 
have obtained the quantum geometric phase of an open 
system undergoing nonunitary evolution. The geomet- 
ric phase is computed in a tomographic manner, i.e. we 
measure the off-diagonal elements of the reduced density 
matrix of the system, from which we extract the decoher- 
ence factor that we use in the definition of the open sys- 
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FIG. 2; (a) Gate sequence for preparing the pseudo-pure state 
|00)5£ from the thermal state pth ~ "fcZs +^hZ£, where jc 
and jH denote the gyromagnetic ratio of C^^ and If^. (b) 
Gate sequence for the quantum simulation of the system and 
measurement of the decoherence factor r(r), which is pro- 
portional to (2(T+) = {X + iY). In both plots the rectan- 
gles denote single-qubit gates, implemented through radio- 
frequency pulses. The rotation angle is shown inside the rect- 
angle, and the direction above. In the experiment we used 
cos 8 — ~ 1/2 and tana — —A/B with a G (0,7r). 

The narrow black rectangles denote the gradient pulses along 
Z-axis. The two filled circles connected by a line denote the 
J— coupling evolution e""^^'S'^£ , where (j) is indicated next to 
the line. 



tem GP. In future work, we will introduce a third (probe) 
spin to perform an independent and direct measurement 
of the GP using traditional interferometry-based tech- 
niques [5] . Our experiments support the observation that 
when the environment is near a second order quantum 
phase transition, the correction to the GP becomes sin- 
gular. For our experiments, we introduced a simplified 
two-level model that captures the essence of the spectral 
behavior of the critical environment: the closing of the 
gap. By adding stochastic fields and further spins, we can 
quantum-simulate more realistic environments and cou- 
plings to the system, so that each initial state is affected 
differently by the bath. Despite the apparent simphcity 
of our experiment, we believe that the techniques we de- 
veloped are quite general and applicable to more complex 
quantum simulations, and to related approaches such as 
bath engineering |20j — designing an environment so that 
it induces a system to relax and decohere to interesting 
quantum many-body pure states. 
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Appendix A: Supplementary Material 

In these supplementary notes we show how the open 
system geometric phase (GP) behaves in the hmit of weak 
couphng to the environment. After analyzing the struc- 
ture of the dominant terms, we will obtain analytically 
closed formulas for the case of an Ising spin chain envi- 
ronment, and compare to exact numerical results. The 
analytical results will show explicitly the singularity of 
the GP when the environment is at the critical point of 
a quantum phase transition. 



a. Small coupling expansion of the geometric phase 
The GP for an open system [Eq. (3) in the main text] 



sin(0+/2) 



|r(t)|sin(0) 



^|r(OPsin2(0)+4(e+-sin2(0/2))^ 



= ,(A3) 



and the only relevant eigenvalue of the reduced density 
matrix of the system is 



£+ = 



\(^ + \/cos2(6l) + |r|2sin2(6l)^ . (A4) 



IS 



2 



tan 

where 
cos(0+/2) 



sin(y3(r) sin(^) sin(|) 



cosi^(t) sin(^) sin(|) -I- cos(^) cos(|) 



(Al) 



We want to obtain a more tractable expression of the GP 
in the limit of small coupling between the system and the 
environment. For this, we expand the decoherence factor 
r{t) = |r(i)|e"''''' in powers of the system-environment 
coupling strength 5, 

|r(i)|2 = l-i?(2)(t)52-i?(3)(i)'5' + 0(<54) (A5) 

The zeroth order in tp can be assimilated as an overall 
phase in the environment, and the first order in |r(i)|2 is 
zero for environments with a finite spectral band width, 
as with spin environments [5 . At this level of approxima- 



2(e+ -sin2(0/2)) 



(A2) 



tion, the arc-tangent term in Eq. (Al) can be expanded 
as 



while the term with the integral is 



V{3) (r) - - cos 6* cos (t) I i?(2) (t) + ^ 



(A6) 



TT (1 — cos6') — sin ^f{i){T)S — 



n 



cos sm 



dt i?(2)(i) 



sm^^-^^3){T) + loose sm^e(nl^ dtR^3){t)-l^ dt R(2){t)^^{t)^ 
where we have assumed that ip(^i^{0) — 1^(2) (0) — (^(3)(0) = 0. Adding up the two contributions results in 



(A7) 



$ ~ Trfl — cos ( 



s^]n f i?(2)(t) dt 

4 Jo 



+ (3R(2){r)^ii){r) + ^l^{T)+6n [ R^3j{t) dt ~ 6 [ R^^)(t)^{t) dt 



(A8) 



h. GP from an Ising spin chain environment 



with a homogeneous transverse field, with Hamiltonian 



Let us consider as the environment a paradigmatic ex- 
ample of quantum criticality: the Ising spin chain model 



HsiX) = -J Z^Z,,+, + XX,}j 



(A9) 



6 



where N is the number of spins in the chain, J the spin- 
spin coupUng, A the dimensionless strength of the exter- 
nal field, and X„ and Z„ are the Pauli matrices of the 
n-th spin. The quantum critical point is at Ac = 1 [3]. If 
the system spin couples homogeneously to the Ising chain 
with strength 6 (i.e. Hgg = S'^^^Zn) this model can 
be solved analytically using a standard Jordan- Wigner 
transformation [1, into free fermionic modes. Notice that 
the requirement of homogeneous coupling is only for sim- 
plicity and can be relaxed in general [4j. The decoher- 
ence factor induced by this environment — with the chain 
initially in the ground state of — decomposes into a 
product of factors, each coming independently from a 
different bath mode P, 



r(t) = Y[Rk{t)e 



i{ipk{t)-ekt) 



(AlO) 



fc>0 



where £k = 2-\/l + A^ — 2Acos(fc), with k = 

This particular product form 



_i_i_^2L -4- H -4- 



. 3 27r 



JV-1 27T 



2 N 

stems from the fact that the environment modes are non- 
interacting, each contributing and independent factor 



fk{t) 




(All) 



J. I cos £kt + i sin ikt cos 2ak 
cos ikt — i sin ikt cos 2ak 



where ik = 2^1 + (A + - 2(A + S) cos(fc), 2ak = 
[^.(A + 5)- ek{X)i and tan(0,) = 

In order to use the approximate expression Eq. ( A8 1 , 
we now expand each factor R^ (t) and (pk (t) in powers of 
the coupling strength (5, 

Rl = 1 - RkX2)it)6^ - Rk.i3)it)S^ + 0{5^) 
Mt) = ekt + 'fk,{i){t)5 + 0{6^), (A12) 

and insert them in Eq. ( |A10[ ), 

= llRlo,ll[l-Rk,i2){t)S'-Rk,i3){t)S'] 



k>0 k>0 

l-52^i?,,(2)(t)-<53^i?fc,(3)(i) 
fc>0 fc>0 
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Rk.(2)it) dk ^ S"" 
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Jo 
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(A13) 



where in the last operation we approximate sums over k 
with an integral — which is a good approximation in the 
thermodynamic limit N ~¥ oo — , and with 



Rk,{2) 



16 



sin^ ksiT?{ekt) 




2 -10 -1 -0.1 
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FIG. 3: Correction 5<1> = $ — $0 to the geometric phase of 
a system spin in presence of an Ising chain environment (cir- 
cles). In panels (a) through (d) we show the correction as a 
function of the strength A of the transverse field of the envi- 
ronment chain. The values of the self-energies of the system 
are Q, = IJ, 2 J, 5J, and lOJ, for panels (a), (b), (c), and (d) 
respectively. Here J is the interaction strength between spins 
in the environment. In solid line is the third order approxi- 
mation, and in dashed line the second order one. The phase 
correction is shown normalized by the length of the environ- 
ment chain N and the strength of the coupling to the envi- 
ronment squared, 5^. In all plots A'^ — 100 and 5 = 5 10~^J. 
In panel (e) we show the correction 5$ at the critical point 
of the environment (A = 1), as a function of the self-energy Q 
of the system. Notice that il is inversely proportional to the 
contact time with the environment, r — 2n/Q. The dotted 
horizontal lines indicate the energies that correspond to the 
left panels. 



Rk,i3)it) 



'/ffc,(l)(0 = 



128(cos fc — A) sin^ k sm{ekt) 

^k 

X [sin(£fci) - ektcos{ekt)] 
A — cos k 
£k 



(AM) 



With these coefficients, the time integrals in Eq.(A8) 



can be done analytically using commercial software like 
Mathematica, which gives us 



— cos u sm 



^^2(A) + -(3T/2(A)Gi(A) 



(A15) 
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2^ Jo 8^4 

[iTTEk (2 + cos(2efer)) - 3n sin(2efcT)] , 
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4A 
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(A16) 



where E{x) and K(x) are the complete eUiptic integral 
and the complete elliptic integral of the first kind respec- 
tively. As we can see, the GP of the system spin must 
be singular at the critical point of the bath A = Ac = 1 
because K{x) has a singularity at a; = 1. 

We computed the GP of the system-spin exactly for en- 
vironment chains of up to 100 spins. Longer chains can 



become computationally unstable and only add fine de- 
tails to the singularity around the critical point. We show 
in Fig. ([3]) the correction to the GP induced by the cou- 
pling to the environment as a function of the transverse 
field of the environment, and for different cycle periods 
of the system. Notice in the figure that the second or- 
der approximation to the exact formula performs poorly 
compared to the third order one [Eq. ( |A8| ], especially 
for large periods r where the environment acts for more 
time. As is to be expected, the duration of contact with 
the environment also affects strongly the magnitude of 
the correction to the GP. At the point where the environ- 
ment is critical, we observe that the GP becomes singular 
in the thermodynamical limit — in contrast to the dis- 
continuity of the GP observed for a first order transition 
of the environment [2] . We can see the singularity in the 
thermodynamical limit appear already in the analytical 
approximate expressions obtained from Eq. ( A8 1 . 
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